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ABSTRACT 


A method for identifying and - after identification - 
reducing linear time invariant systems on the basis of input - 
output records of the system is reviewed, applied to physical 
devices and extended to handle the case where noise is present. 
The method is implemented using an analog to digital converter 
in order to have the input and output records in a digitized 
form and thus to be able to use a digital computer program 
composed of a numerical integration subroutine and another 
subroutine to solve overdetermined sets of linear algebraic 
equations. Several practical examples are presented to 
demonstrate the accuracy and present capabilities of the 


procedure. 
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I. INTRODUCTION 


A. THE IDENTIFICATION PROBLEM 

All the modern techniques of analysis and design of 
any kind of system are based on the previous knowledge of 
a mathematical model of the system to be studied. This 
principle is also applicable to control systems. The form 
of the model depends on the methods of analysis and design 
to be employed as well as on the physical characteristics 
of the system and since most of the theory on control 
systems is based either on the state space or transform 
representation of systems, the majority of the mathematical 
models will be a set of state equations or a transfer 
function. 

After deciding what basic form the mathematical model 
should take, the problem of determining the numerical values 
of the parameters arises. These values can often be deter- 
mined from the laws of physics and the data supplied by the 
manufacturer or obtained through testing. This is not always 
the case, however; sometimes the laws of physics become 
mathematically intractable or are very difficult to apply and 
quite often the values of certain key parameters are not 
available or have changed with time. It is in these cases 
when the identification problem arises. This problem consists 
in determining the system model given the input, output and 


initial conditions. 





For the purpose of this work, the identification problem 
can be stated as follows: 

Given a system considered a black box, obtain a record 
of the input and output of the system over some finite 
period of time - 400 msecs. in this particular case - and 
from it find a mathematical model - the transfer function 
relating output and input in this case - and the values of 
the parameters in such a way that the model will accurately 
describe the system. 

It should be kept in mind, however, that the problem of 
identifying a system solely on the basis of input and output 
quantities is very rare in engineering, because in most of the 
cases the engineer has some idea whether the system is linear 
or non-linear, time varying or time invariant, what is the 
order of the system itself and perhaps what are the time 
constants involved and in the presence of noise the statis- 
tical characteristics of it. For these reasons rather than 
having a BLACK BOX, what the engineers have is what could be 
called - a GREY BOX -, and as will be shown later, the 
identification techniques depend very heavily on knowledge 
of the characteristics mentioned above and in this particular 
case THE ORDER OF THE SYSTEM. 

Although the control systems literature on system identi- 
fication techniques is quite vast, there are no known identi- 
fication schemes capable of handling all identification 
problems effectively. Some of the proposed identification 


methods are mentioned in the following: 
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Nieman, Fisher and Seborg [1] have pusblished a paper 
summarising most of the common approaches to the problem of 
identifying lumped parameter systems. Eykhoff [2] has a good 
discussion of the industrial application of several methods. 
Since the response to an impulse gives directly the transfer 
Fiction of linear systems, obtaining this response is a 
common scheme used in linear identification. Mishkin and 
Haddad [3] have developed a technique for finding the impulse 
response based on samples of the system output and its deri- 
vatives. Levin [4] and Kerr and Surber [5] have developed a 
scheme for estimation of the impulse response from the 
samples of a noisy input and output. Turin [6] and Lichtenberger 
[7, 8] have used a matched filter to obtain an identification. 
Hill and McMurtry [9] have suggested the use of white noise, 
a binary test signal and crosscorrelation techniques to 
attack the problem. Cooper and Lindenlaub [10] have investi- 
gated the use of matched filter and crosscorrelation. 

Another common approach to the problem is to determine 
the coefficients of the differential or difference equation 
which describes the system. Kumar and Sridhar [11] have 
employed the method of quasilinearization with some success. 
Nagumo and Noda [12] have developed a learning approach to 
the problem. Bass [13] uses modulating functions and the 
method works well in the presence of noise. Astron and 
Bohlin [14] use a statistical optimal method of determining 
differential equation parameters known as the "maximum 


likelihood method". Peterka and Smuk [15, 16] have developed 





a similar method which is not optimal but is considerably 
simpler computationally and it is known as the "prior 
knowledge fitting method". Ho and Kalman [17] have proposed 
an algorithm for determining state variable models of sampled 
data systems which works well in the presence of noise, and 
has Te extended to continuous systems by Eldem [18]. 
Identification methods vary widely with respect to how 
much must be known about the system before the method can be 
applied. In general, the less known about a given system, 
the more difficult it is to find a method capable of accom- 


plishing the identification. 


B. OBJECTIVES OF THIS PAPER 

This paper will present a research, handling actual 
hardware, of an identification technique originally suggested 
by Diamesis [19] and studied using simulation techniques by 
W. R. Hansell [20]; also an attempt at using the technique 
in the presence of noise is made. The method TT a 
knowledge of the system input and output over some finite 
time interval and a rough estimate of the order of the system 
is desired. Although, theoretically the system input used 
need not be restricted to a class of testing functions, that 
is, it can be completely arbitrary, experimental results 
showed that inputs such as step, ramp or in general any 
periodic wave, sine, square, etc., should be avoided. 

Unlike some identification schemes which require the 


calculation of derivatives of the input and output, the 








technique to be presented requires only integrals of the 
input and output, and the advantages of numerical integration 
over differentiation are well known. 

In Chapter II, the theoretical development of the 
identification technique is given. A method for identifying 
the Sc conditions of the unknown system is also presented, 
but in this case it is necessary to know the value of the 
derivatives of the input at time zero. Since the initial 
conditions do not modify in any sense the mathematical model, 
their knowledge is not important. 

Chapter III presents the method used to check the 
developed theory in Chapter II. Several examples are given 
to demonstrate the capabilities of the method. 

In. Chapter IV, recommendation and conclusions are given 


and following this chapter a complete listing of all computer 


programs are presented. 
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II. IDENTIFICATION BY MULTIPLE INTEGRATIONS 


A. GENERAL APPROACH 

The development that follows is similar to the ones 
given by Diamesis [19] and W. R. Hansell [20]. The develop- 
ment given by Diamesis is restricted to the case where all 
initial conditions are zero and the parameters of the model 
are the unknowns of a uniquely determined set of linear 
algebraic equations. Hansell extended the development to 
the case where the initial conditions need not be zero and 
he makes use of an overdetermined set of linear algebriac 
equations with the model parameters and some coefficients 
which have information about the initial conditions, as 
unknowns. 

The overdetermined set of equations are solved using 
the method of least squares. The method is extended to 
handle noisy inputs and outputs. 

Any single input, single output, lumped parameter, linear 
time invariant system can be described by a linear ordinary 
differential equation with constant coefficients. The 
general form of this equation is given in equation (1). A 
set of initial conditions is also necessary to completely 
define the system. 


y™ (t) та уЗ (е) аа а а т a y (t) - 


n-1 
(1) 


= Ь 0 (е) Eccc et b U(t) where, in a general case 


n>m and 


10 


(3 апа 02) (=) represent the a and й derivatives 
of the output and input, respectively. 

The necessary initial conditions are: 

The value at time zero, of the en and its n-1 first 
derivatives and the value at time zero, of the input and its 
m-1 first derivatives. 

The notation to be used is: 

У (о), yt! (o) M rM У’ (о), y Em qoe to denote 
the, n-1, І.С. of the output and 

u(o), а (о) И е ж TE “2 (6), es to denote the, 
mel, IC. On the input. 

The identification problem to be treated here consists of 
determining the coefficients а. апа Do of equation (1), on 
the basis of input and output records taken over some arbi- 
trary time interval, for every Palmar values chosen forn 
and m. The initial conditions will be assumed to be unknown 
and also it is assumed that the noise present is zero mean 
and that the standard deviation is small compared with the 
R.M.S. of the input and output, for the technique to work 
properly. 

Taking the Laplace transform of (1), yields equation (2) 
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where the d. coefficients account foxr the contribution of 
the initial conditions. 


Dividing equation (2) by Aue results in equation (3) 
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equation (4). 
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And since the system is assumed to be time invariant, 


nothing has been lost by setting the lower limit on the 


integrals equal to zero. 


Rearranging terms in equation (4), equation (5) is 


obtained. 
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Since records of the input and output are assumed to be 
known, a linear algebraic equation with the system parameters 


а;, da and the terms that contain information about the 


initial conditions as unknowns, can be formulated by performing 


the multiple integrations, from zero to some time t or from 


EN 


any time to to a final time t after assuming some values 


i 


for n, and m. A set of 2n+m+l equations can be obtained by 


№3 





letting t. take on 2n4*m*l different values. Assuming that 


f 
the equations are linearly independent, it will now be 
possible to solve for the 2n+m+l unknowns; n+m+l, а., z 


ce frficrents and the n а, terms which contain information 
about the initial conditions. It will be shown later how 
Е 2 coefficients can be used to determine the initial 
conditions. 

Consider now the case when noise is present at the input, 


at the output, or at both. 


Under these conditions equation (5) becomes: 
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dtf = - y(t)+v(t) dt 


where v(t) and r(t) represent the noise at the output and 
input respectively. 

If the coefficients in equation (6), that is, the 
multiple integrals of the input, output, and time, are going 
to have the same values as the corresponding ones in equation 
(5), in order to have the same solution for the asi, Dae and 
Ч, in other words, the same identification, the presence of 
noise should not affect the value of the integrals. 

Consider a general term: 

Я E 
COL [ree] at? (7) 
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Equation (7) can be written: 
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The first integral of equation (8). 


te t 
eee (p) esė y (t)atP (9) 
O O 
is the true value of the coefficient of sal in equation 
(5) ата 
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^T" & Ў... v (t)atP (10) 


О О 
the undesired contribution of the noise. Thus, what has to 
be done is to find when the value given by equation (10) is 
negligible with respect to the value of equation (9). 
Assuming zero mean noise, the expected value of 
equation (10) is: 


te t 
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Since zero mean noise was assumed and the expected value 
of the integral is equal to the integral of the mean value. 
But the fact that the mean value of equation (10) is zero is 
not enough to say that the influence of the value of equation 
(10) in (7) is negligible. Besides the mean of equation (10), 
its standard deviation or R.M.S. value should be found in 
order to be able to compare it with equation (9), and if the 
R.M.S. value of equation (7) is very small compared with the 
value of equation (9), it can be said that the influence of 
the noise will not significantly affect the identification. 

The variance of equation (10) considering, р - 2; Tor 


simplicity, is given by equation (12) 
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where ү апа £ are dummy variables. Equation (12) can be 


YA [| У (1 ) Үү ds; / pe = 
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but E Ku is recognized to be the auto-correlation and 





written 


can be written as R(t), where T Eos В 
The R.M.S. value of equation (10) can now be expressed 


as follows: 


ae ин 1 
R(t) dy, dyj,dé,d&, 2 С 
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1 
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E (о) ауаё ; В, (т) < В, (0) (13а) 
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or going back to the initial notation 


te t 
E VR, (o) at* (14) 
O O 


If now (l4) is compared with (9), with p = 2, it is observed 
that both integrals have the same form. Therefore, if 

y(t) >> R. (o) = R.M.S. value of v(t), the noise is not 
going to have much influence in the identification. The 
same conclusion can be reached for any value of p. 

As was said before, the set of 2n+m+l linearly indepen- 


dent algebraic equations obtained from equation (5) letting 


t, take 2n+m+1 different values, would solve the problem of 


T 


identifying the parameters а., ©, апа а but since for 


some kind of inputs and some t the 2n+m+l equations may 


Е’ 
not be linearly independent, thus instead of forming just 
2n+m+1 equations, a set is formed which has a number of 
equations substantially bigger than 2n+m+1 and this over- 
determined set is solved using the method of least squares. 

This technique also allows some flexibility in the 
Selection of the initial values of n and m as will be shown 
later. 

If in equation (2) the initial conditions are set equal 


to zero, after a simple manipulation, the transfer function 


of the system can be obtained and is given in equation (15). 





Y(s bs +... + b.S + b 
= M I O (15) 
n n-1 : 
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It should be noticed that the coefficients а, апа Б. 
of equation (10) are exactly the same as the а; апа Б. of 
equation (5). Therefore, as soon as the set of overdetermined 
algebraic equations is solved for a, and re the transfer 


function is determined. 


Equation (15) is usually written in the form (16). 


(16) 








Where K = ба is the so called constant gain and Cg = 29 for 


m 
Cas g <m. 


If the roots of the numerator and denominator polynomials 
are found, the more familiar form of the transfer function 


can be written; equation (17). 


U (s) к (5 + р;) 
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Where -р; апа БИ are the poles and zeros respectively. 


Equation (16) can also very easily be written in state 


form as 
i 
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But mhowsare » and m determined? Те мас зала that im 
order to apply the method efficiently a rough idea of the 
order of the system should be available, although it is not 
absolutely necessary; that rough idea is used to pick n' and 
m' as the first approximation for n and m, taking into account 
that n'>n and m'>m to be sure that neither poles nor zeros 
are missing Since m iseusually smaller thamn, m' should 
not be greater than n'-l. 

Consider then that n and m are the true numbers of poles 
and zeros and that n' and m' are assumed. Under these 
conditions one would expect to obtain n'-n, or m'-m - what- 
ever number is smaller - poles equal to zeros, what is 
equivalent to saying thai they Bec not present in the system, 
what is actually the case; the (п' - п) - (п' - m) poles or 
(n'- n) - (n' - n) zeros in excess which have not been 
cancelled will appear to be either in He right half plane 
or with values too big compared with the actual values of the 
poles and zeros in the system. As a result after a first 
guessing for n and m, a clear indication of the actual order 
of the system and of the number of zeros, should show up. 


Examples in Chapter III will demonstrate this point. 
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B. IDENTIFICATION OF THE INITIAL CONDITIONS 

Although the knowledge of the initial conditions is 
not important by itself, it might be if one desired to 
compare the system model with the actual system by exciting 
the system model with the same input that was used in the 
identification. By comparing the output of the model with 
the output of the system a rough idea of the accuracy of 
the model can be obtained. However, in order for this to be 
possible, all the initial conditions of the input and output 
should be known. 

When the overdetermined set of linear algebraic equations 
is solved, besides the values of а. апа BER the d, initial 


k 
condition terms are also found. These d. terms are related 
with the initial conditions by equation (18). This equation 
can eaSily be found by specifically writing in the contribu- 


tion of the individual initial conditions when the Leplace 


Transform of equation (1) was taken, to obtain equation (2). 


n-1/ Y (0) 
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(18) 
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Equation (18) requires the knowledge of m-l derivatives 
of the input and numerical techniques may be necessary to 
find them, however, this difficulty may be avoided in many 
cases by choosing the beginning of the input-output record 
at a moment where the input is constant and thus all the 


derivatives would be zero. 


С. SPECIAL CASES WHERE SIMPLIFICATION COULD BE POSSIBLE 

Theoreticaily the method should a the -system 
independently of the number of unknowns,but practically 
this is not the case due to the natural limitations in the 
accuracy of the numerical integration methods and in the 
measure of the input and output records, thus it would be 
desirable to reduce the number of unknowns as much as 
possible, because as a general rule the fewer the unknowns, 
the more accurate its identification. 

A great reduction in the number of unknowns occurs if 
the record of the input and output begins when all the 


initial conditions are zero. Under these conditions, all 
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the terms d, of equation (5) are zero and the number of 


K 
unknowns is reduced from 2n+m+1 to n+m+1 since it was 
assumed n>m. 

Any additional information could be used to simplify 
the problem, for instance if the value of the steady state 


gain constant were known the number of unknowns could be 


reduced by one. 
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III. IMPLEMENTATION 


a 


A. NUMERICAL METHODS 
As was mentioned before, after having the record of the 

input and output of the system, which can be obtained as 

described in Appendix I, it is necessary to perform the 

multiple integrations indicated in equation (5), for 

different times te; in order to form the set of overdetermined 

linear algebraic equations whose solution gives the values 


к” 


are of interest in this work. The а 


of the parameters ai, ur and d although only the а; апа БЕ 


кұ, even available are 

not used. _ 
The multiple integrations can ђе handled for any numerical 

method subroutine and although there are quite a few highly 

sophisticated numerical integration procedures available, 

trapezoidal integration will be used because most of the 

more complex integration techniques perform poorly when the 

function being integrated is discontinuous. Since the 

inputs used here are step-wise waves, complex integration 


techniques are not desirable. On the other hand, the slight 


improvement in accuracy offered by more advanced methods is 


not, in general, enough to justify the tremendous increase 
in computational load associated with their use. 

The solution of the set of overdetermined linear equations 
is a classical problem in several fields of mathematics and 


engineering. Most of the classical techniques for solving 
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m = 





such problem are not practical. They tend to magnify 

the errors introduced by the finite precision of the computer 
to the point where the solution is meaningless. Fortunately, 
several modern methods are available that display more 
acceptable behaviour. The method used here was developed 

by Golub [22]. The basic approach is to triangularize the 
coefficient matrix by performing a Choleski decomposition. 
The decomposition is accomplished by applying repeated 
Householder transformations [23]. Once the coefficient 
matrix has been triangularized the unknowns can be obtained 
by back substitution. The method is quite stable numerically 


and is capable of handling ill-conditioned coefficient 


matrices. 


B. CHARACTERISTICS OF GOOD INPUT-OUTPUT RECORDS 

The accuracy with which a system can be identified is 
strongly dependent on the input-output record used in the 
identification. Since parameters are identified on the 
basis of their effect on the output, it will be impossible 
to identify a parameter unless its effect is measurable. If 
the input to a system has a frequency spectrum that is more 
or less uniform over the frequency range of interest, the 
identification will probably be very good. If the frequency 
Spectrum of the input is confined to a narrow band of fre- 
quencies, the identification will probably be very bad. It 
is well known that signals with sharp discontinuities have 


a broader bandwidth than slowly varying continuous signals. 
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For this reason, the step-wise wave was chosen as the input 
in this research, after having demonstrated experimentally 
the preceeding observations seeing that with a step input 

it was never possible to obtain a good identification even 
for a second order system. Periodic functions should also 

be avoided at the input, because very probably many of the 
algebraic equations are going to be linearly dependent. This 
fact was also proved experimentally; a periodic square wave 
was tried as input and the identification was also bad when 

a step-wise input gave very good results. 

Another desired characteristics of the input, experimentally 
ваем, іс ав it should Бе adjusted in such а way that the 
output, if it is strictly increasing, should at some moment 
in the sampling interval have a negative second derivative; 
in other words, the output is desired to have an inflexion 
point within the period in which the record is taken. The 
reason for this is thought to lie in the subroutine DLLSQ 
which solves the overdetermined set of algebraic equations. 
What DLLSQ essentially does is look for a fitting output 
function that minimizes the least Sonate error between the 
generated and the measured and if the second and higher 
order derivatives of the output are positive at the end of 
the sampling interval, the function generated by the sub- 
routine may in some cases correspond to an unstable situation. 
Examples in Chapter IV will also ee UNITS point., 

In relation to the sampling rate at which the input and output 


should be digitized, it can be said that low sampling rates 


26 





may introduce inaccuracies in the. numerical integrations 
and as a result poor identification, especially in the 
location of high frequency poles and zeros. Maximum 
sampling rate available or ten to one hundred times shorter 
than the shortest system time constant is recommended. The 
output of the system is also recommended to have maximum 
value of the order of the maximum value of the input, thus 
the numbers with which the computer deals are not very 
different in magnitude. If the output has not the same order 
of magnitude of the input, it can as a rule be amplified 
accordingly, taking this amplification factor into account 


when the time comes to see what the actual identification is. 


С. PROGRAMS DESCRIPTION 

Four different computer programs were used in this 
Merch. A eine of them is given after Chapter IV. The 
first one is the A/D program. This program is used for high 
speed A/D (analog to digital) conversion. Its output 15 a 
seven track magnetic tape that contains the digitized input 
and output record of the system. This program consists of 
a main driver program and five subroutines. 

The source language is FORTRAN for the main program and 
METASYMBOL for the subroutines. 

The A/D program provides a means of sampling and digitizing 
analog information at a rate (up to 10KHz) determined by an 
interrupt signal originated at the analog computer logic 


patchboard (Ci 5000). (See Appendix I). The digital data 
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are stored sequentially in one of two buffers provided by 
the main program. When one buffer is full, the storage 
sequence is continued at the top of the other buffer. An 
indicator notifies the main program of this switching and 
the main program processes the data in the full buffer 
ыа them in the magnetic tape while the other buffer 
is being filled. 

The parameters to be entered in the main program are: 

NSAMP = Number of words in: one block/NCHAN. 

NCHAN = Number of channels to be digitized. 

NREC = Maximum number of records to be digitized in 
one file. 

ITAPE = Number of tape drive unit to be used. 

NDEL = Delay between digital and analog signal in 
option 7. The options of this program are: 

1 = Desired to input new parameters. 

2 = Start digitizing the analog input signals. 

3 = Write an END OF FILE on the magnetic tape. 

4 = Rewind the tape to the loading point (beginning). 

5 = Skip the following number of END OF FILE marks. 

6 = Print the following amount of digitized data. 

7 = Actuate the A/D program for the next single block 
of data encountered and feed the signal to the recorder. 

The five subroutines are: ADSTART, ADSTOP, ADFAST, 
MTRDY and MTOUT. 

Subroutine ADSTART initializes the A/D input operations. 


The calling sequence is: 
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CALL ADSTART (NCH, BUF, NEWBUF, NSAM, RECNUM, NEXLOC), 
where 


NCH = Number of A/D channels to sample. 


BUF A pointer to the first buffer. 

NEWBUF = A pointer to the second buffer. 

NSAM = Number of samples per buffer. 

RECNUM = A record counter which is bumped by ADFAST 
each time a buffer is filled. 

NEXLOC = The return address for ADFAST to be used each 
time a buffer is filled. 

The first four arguments are used to set up pointers 
and counters in ADFAST.  RECNUM provides a means for ADFAST 
to keep the main program aware of the record number of the 
buffer currently being filled.  NEXLOC is an alternate 
return location for ADFAST.  ADSTART also sets up the 
interrupt controller branch instructions and starts the A/D 
hardware.  ADSTOP stops the A/D operations and restores the 
interrupt locations to their normal instructions. 

ADFAST controls the buffer location pointer which directs 
the A/D data into core. 

MTRDY initializes the magnetic tape handler. It has 
the following arguments: 


UNIT - Magnetic tape unit number. 


BUFl - First word address of first buffer. 

BUF2 - First word address of second buffer. 

NW = Number of words to be written per record. 
IR = Tape operation status word return. 
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1 if tape I/O channel is busy. 
2 if operation completed. 

3 if tape unit is not ready. 

4 if tape writes error. 

MTOUT controls the stream of A/D data from the buffers 
into the tape. It is called by the main program each time 
a record of data is ready for recording. The argument, NB, 
indicates which buffer is just eee Ama is to be output. 
After setting the status indicator (IR) to 1, to indicate 
an operation not complete condition, the routine tests the 
tape and channel status. If either is busy, IR is bumped to 
(not er condition) and the routine exits back to the 
calling program. If everything is all right, NB is checked 
and the I/O instructions for writing the filled buffer on 
tape are executed. If the operation is completed without 
error, IR. is set to 2 (operation completed) and to 4 (write 
error) if an error occurred, and execution is returned to 
main program. 

The seven track magnetic tape generated by the A/D 
program recording the digitized input and output can not be 
directly read by the IBM 360 computer system, and since this 
system is used to perform the integrations of equation (5), 
a seven to nine track tape conversion should be done prior 
to the numerical integration. This conversion is done by 
the IBM 360 using the CONVERT program, which is described 
in Ref. 21. This program has several parameters. One of 
them is BLKSIZE which is related with the NSAMP and NCHAN 
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parameters used in the A/D program.  BLKSIZE can not be 
greater than 32.700 and for N€HAN - 2 limits the maximum 
number of samples that can be taken for the input and out- 
put record to about 4000. 

The third program is the main one and is a direct 
implementation of the procedure developed in Chapter II. 
This program has several important parameters that have to 
-be defined at the very beginning. NP and NZ are rough 
Enemuaebesvossehewmumber Ol Poles and zeros" in the System to 
be identified. NP and NZ stand for n' and m' in the program, 
N' and m' were defined in Chapter II.  KPTMAX is the number 
of sample points available in the input-output ње m The 
input amplitude is R and the output amplitude is C. Because 
the integration of the time has to be done, the program 
generates the time corresponding to each sample for what 
the sampling rate has to be known and is represented by the 
parameter DT in the program. IPTS is the number of sample 
points that will separate successive linear equations. 
Therefore, the total number of linear algebraic equations 
formed is equal to KPTMAX/IPTS. The program reads: in the 
nine track tape generated by the CONVERT program, performs 
the integrations of equation (5) and forms the overdetermined 
set of linear algebraic equations. When the set is formed, 
subroutine DLLSQ is called. This subroutine is an implemen- 
tation of the Golub algorithm for solving overdetermined 
sets of linear equations.  Subroutine DLLSQ returns the 


values as, Ds and Ч. Тре а, апа De coefficients are fed 
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into the RTPLSB subroutine, which is a polynomial root 
finder, in order to obtain the poles and zeros of the system. 
The output of the program are the results of the 
identification given in both transfer function and state 
variable form. 
The fourth program was written to compare the outputs 
of the actual system and of the ones obtained from the 
identification, when all of them are excited with a step. 
This program solves the differential equations that describe 
the system. The differentials equations should be arranged 
ме сасе variable form. “The output of the program is a 
drawing of the response of the different systems to be 


compared in the same plot. 


D. EXAMPLES 

The following examples show the possibilities of the 
method, how the order of the system can in general be 
determined, if it is not already known, from a trial identi- 
fication run using an estimated order greater than the 
actual order of the system. They also show how the accuracy 
of an identification depends on the input wave form, on the 
R.M.S. value of the noise present, and as expected, does 
not depend on the initial conditions. They also indicate 
how the method can be used to reduce or obtain a lower order 
mathematical model. In all the examples, the MAIN PROGRAM 
PARAMETERS were: 


NUMBER OF DATA POINTS = 4000 
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NUMBER OF EQUATIONS 40 


100 


|| 


DATA POINTS PER EQUATION 

In examples 1 through 6, the input function was a step- 
wise increasing wave form um form, see Fig. 1, 
Example 1), with a maximum value of about 85 volts. The 
amplitude and the number of steps were changed in each 
example. In Example 3 and 4,the maximum amplitude of the 
input was about 9 volts. In examples 6 through 10, the 
input function was also a step-wise wave form, but instead 
of being strictly increasing, it had one part strictly 
increasing and the final part strictly decreasing. See 
Fig. 4, Example 5. | 

The systems investigated in Examples 1, 2, 5, 6, 7, 8, 
9, 10 and 11 were set up using R-C and R-C-L circuits 
isolated by operational amplifiers which also provided the 
desired amplification factor. In Examples 3 and 4, an 
active system was studied. This system was a D.C. positioning 
Servo. | 

Example 1 was performed as a first test and in ita 
second order system - two real poles - was considered. The 
input was a step-wise increasing function, (see Fig. l), 
with zero initial conditions. The results of the identifica- 
tion were within the 2% for the dominant pole and the gain 
constant and within the 10% for the higher frequency pole. 
This example shows how assuming different values for m' and 
n', the results of the identification gave a clear indication 


of what the actual order of the system was. 
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Example 2 illustrates the identification of a third 
order system. The input function as in Example 1 was a 
step-wise increasing function. The initial conditions were 
also zero. The identification of the low frequency poles 
and the gain constant were also within 3% of the true values; 
however, the pole at -100.0, was identified with an error of 
about 30%, this result was not a surprise since the more 
negative the value of a pole, compared with the other poles 
in the system, the less exact its identification. This 
example shows also the fact that in some cases the trial 
runs tO Lind out the true values ot n' and m', did not indicate 
very clearly what these values were. In cases like this, a 
thorough search should be done around the expected values of 
n and m and different input-output records are recommended 
Be be used. 

This example also shows how the method could be used to 
reduce the order of the system, i.e., when п' = 3 and m' = 
l were assumed, the results correspond practically to a 
second order system, since the pole at -2.22 and the zero at 
-2.20 cancel each other. When n' = 2 and m' = O were 
assumed, practically the same results were obtained. The 
steady state error coefficient of the second order has 
practically the same value of the third order. Essentially, 
the method neglected the effects of the higher frequency 
pore. 

The experiment described in Example 3 was performed 


using a D.C. positioning servo. The transfer function of 
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the servo was first obtained from the asymptotic approxima- 
tion of the closed loop frequency response plot shown in 
Fig. 2. The parameters of the servo were adjusted in such 
a way that two dominant complex roots were present. During 
the frequency response test, it was found that the servo 
had nonlinearities since a jump in the resonance frequency 
was observed and also that reducing the amplitude of the 
input sine wave from twenty volts to one volt peak to peak 
the resonance frequency was shifted from 2.6Hz to 8.1Hz 

(52 rad/sec). 

In taking the input-output record in this example, the 
input was adjusted to be a step-wise increasing function, 
(see Fig. 1), but with a step amplitude of about 0.5 volts. 
The 4 volts step amplitude used in Examples 1l and 2 was too 
big; the D.C. motor ran continuously. 

The results of the identification were almost identical 
to the ones obtained from the frequency response plot. 

In Example 4, the same D.C. servo was treated but now 
the parameters of the servo were adjusted in order to have 
two real poles. Also a frequency response test was performed 
and a transfer function was obtained from the frequency 
response curve asymptotic approximation (see Fig. 3). The 
results of the identification were very close to the ones 
found from the frequency response plot. 

In Example 5, a third order system is again considered, 


but now a zero was added. 
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тће input function and the I.C. were the same as in 
Example 1. The results in this case also showed a poor 
ые івісаеіоп of the higher frequency pole. 

Example 6 shows the effect of the input function in 
the accuracy of the identification. Now a fourth order 
system is studied uSing two different inputs. The first 
input was the same as in Example 1 and is shown together 
meen the output of the system in Fig. 4. With this input, 
the results of the identification were very poor since not 
a Single pole was correctly identified. Several fourth 
order systems were treated with this type of input and the 
identification was bad, giving in some cases poles in the 
right half plane. It was observed, however, that the output 
of the system, at the end of the record had a positive 
second derivative" (see Fig. 4). It was thought that perhaps 
this fact could be the cause of the poor identification 
considering that an unstable system might have the same 
output, for the finite interval of time considered, and this 
unstable system could be the one identified or generated 
by the DLLSQ subroutine. With the preceding in mind, the 
input function was changed and adjusted in such a way that 
the output had the second derivative negative at a certain 
instant within the time interval considered (see Fig. 5). 
The results in this case were very good since even the high 
frequency pole was identified with an error within the 6$ 


of the true value. The I.C. were also zero value. 
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In Example 7, the identification of a system with four 
poles and two zeroes is presented. The input function was 
as the one shown in Fig. 5. The results in this case were 
also very good. 

Example 8 shows how the I.C. does not affect the identi- 
fication, the results being practically the same with I.C. 
equal to zero as with І.С. 0; several different initial 
conditions were considered. 

In Example 9, a third order system was studied, but noise 
was added at the input and output. The results showed the 
great influence of the noise in the identification. It was 
found that in order to have a good identification, the 
maximum R.M.S. value of the nois@ should not be greater than 
about .5 volts. 1.0 volts R.M.S. noise made the identifica- 
tion very bad. The noise added was a gaussian, low frequency 
noise - from a few cycles per second to about 300 cycles per 
second - with.zero mean, 

Examples 10 and 11 show how the method can be used to 
reduce the order of the system. In Example 10, a fifth 
order system was treated. The identification was not good 
for the higher frequency poles. Other fifth order systems 
were studied and in some cases none of the poles was identi- 
fied, which indicates that fourth order systems are perhaps 
the highest order systems that can be identified with the | 
kind of inputs used in this work. | 

In Example 10, the reduction was not,as a direct conse- 


Gence OF the poor Identification, very good as can be seen 


3m 












— | 





D - = 

E a 
a A <=> 
عن‎ ——— o os ам. 
tin axe Ga ^ 
о = — e = — 0 











< 
— — и жы = а 
>> >> и а 
p — ms ж те с 
me > MEE 
=> = зо 
=> „=> Р н 
к” 
= nn 
= ж = 


comparing the responses to a step input of the different 
models, see Fig. 6. However, it is observed that the 


three models generated by the method, сі. s and 2” 


а 
order, give responses fairly close to one another. 

In Example 11, although the identification is considered 
very good, the responses of the different reduced models 
practically coincide during the transient period, but in 
steady state an error of about 30% is observed between the 
step responses of the fourth order identified model and the 
third and second order reduced models (see Fig. 7). However, 
if the gain constant of the reduced models is adjusted in 
order to have the same steady state gain as the fourth order 
model, the responses are practically the same (see Fig. 8). 

It is to be noticed that the systems considered true, 


in the examples, were true within the tolerances of the values 


of the resistors, capacitors, coils and amplifiers used. 
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EXAMPLE I 


SYSTEM TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -10.0 0.0 J 
2 -100.0 | 0.0 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 1000.0 


TESTS FOR EXAMPLE I 


RUN NO. OF POLES NO. OF ZEROES RESULTS 
a 2 | 0 Good 
b Э | 1 See Table of Data 
с 4 2 See Table of Data 
d 2 1 See Table of Data 


INPUT FUNCTION STEP-WISE STRICTLY INCREASING SHOWN IN 
FIG. 1. 


INITIAL CONDITIONS EQUAL ZERO. 
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INPUT 


VOLTS 
80 
о 400 MS ` 
ЛЕППЕН Пр ІП БАДЫРЬ WO. L 
OUTPUT 
VOLTS 
60 
о 


400 MS 


OUTPUT FU: COTON OA PEG 40,1. 


TEGU E L 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
-9.8488868 00 

2 -111.04012 0.0 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 984.74146 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

1 10935 2.0285 

2 120.88900 
B VECTOR 

2 984.74146 
C VECTOR 

1 1.0000000 

REMARKS: 


ASK FOR TWO POLES NO ZEROES. 


GOOD RESULTS. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
+ -97.994740 0.0 J 
2 2.2362954 0.0 J 
3 -9.8646060 0.0 J 
ZEROES REAL IMAGINARY 
1 2.2427425 0.0 J 


GAIN CONSTANT -. 867.90796 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 -2161.7809 
2 725.47414 
3 105.62305 
B VECTOR 
3 867.90796 
C VECTOR 
1 -2.2427425 
2 1.0000000 
REMARKS: 


ASE FOR THREE POLES ONE ZERO. 
RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE SYSTEM 


SINCE ZERO AT 2.24 CANCELS POLE AT 2.23. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 


SYSTEM TRANSFER FUNCTION 


POLES 


iL 


ZEROES 


l 


2 


GAIN CONSTANT 


REAL 
-4.4110354 
-4.4110354 
-99.827986 


-9.8612480 


REAL 
7.090656 74 


-7.0865674 


= 727.2736 


8 


IMAGINARY 
-54.484437 
54.484437 
0.0 


00 


IMAGINARY 
-59.922949 


59.922949 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 


1 


2 


3 


4 


B VECTOR 


4 


C VECTOR 


1 


2 


2941483.4 
336437. 5 
4940.1259 


118.51130 


72 7.275368 


3690. 9793 


142173135 


1.0000000 
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REMARKS: 
ASK FOR 4 POLES TWO ZEROES. 
RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 


SYSTEM, SINCE COMPLEZ ZEROES CANCEL COMPLEX POLES. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (d) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -9.8592476 0.0 J 
2 -97.619069 0.0 . 2 
ZEROES REAL IMAGINARY 
1 -847.73348 0.0 Ј 


- 


GAIN CONSTANT = 1.0221453 


SYSTEM STATE VARIABLES (PHASE FORI) 


A VECTOR 
1 962.47832 
2 -107.47832' 
B VECTOR 
2 1.0221463 
C VECTOR 
Т 847.73348 
2 1.0000000 
REMARKS: 


ASK FOR TWO POLES ONE ZERO. 
RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE SYSTEM. 


SINCE THE ZERO АТ - 847.7 САМ BE NEGLECTED. 
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ЕХАМРЉЕ 2 


Seen TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL ` IMAGINARY 
J -1.000 0.0 
. 2 | -10.00 | 0.0 
3 -100.0 0.0 
ZEROES 
GAIN CONSTANT - 10000.0 
TESTS FOR EXAMPLE 2 
RUN NO. OF POLES | NO. OF ZEROES RESULTS 
a 3 0 Good 
b 4 1 See Table of 
С 4 0 See Table of 
а 4 2 See Table of 
e 3 1 See Table of 
E 3 1 See Table of 


РРБ RÜNCTTON STBP-WISE STRICTLY INCREASED SIMILAR 
TO THE ONE SHOWN IN FIG. 1. 


INITIAL CONDITIONS EQUAL ZERO. 


Data 


Data 


Data 


Data 


Data 





IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTEM TRANSFER FUNCTION 


POLES 
1 
2 


3 


ZEROES 


GAIN CONSTANT 


REAL 


-1. 1072821 
-10.224865 


“ТӘ ФЕН 27 


REAL 


11153.483 


SYSTEM STATE VARIABLES (PHASE 


A VECTOR 
1 
2 
3 

B VECTOR 
3 


G VECTOR 


REMARKS: 


1484.3522 
1429720893 


142.44357 


11153.488 


1.0000000 


IMAGINARY 
970 Ј 
0.0 J 
0.0 J 
IMAGINARY 
FORM) 


ASK FOR THREE POLES NO ZEROES. 


GOOD RESULTS FOR LOW FREQUENCY POLES, NOT SO GOOD FOR 


HIGH EREQUENCY POLE. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
it =1, 2075351 0.0 J 
2 -9.9863699 0.0 J 
3 =1. 76 OOS И, Ј 
ZEROES REAL IMAGINARY 
1 157. 18987 0.0 : J 
GAIN CONSTANT = -1919.7240 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
р 41903.747 
2 40128.130 
3 3645.1011 
4 14.692118 
B VECTOR 
4 =1919,. 1240 
C VECTOR 
d: -1579918907 
2 1.0000000 
REMARKS: 


ASK FOR FOUR POLES ONE ZERO, 
RESULTES ТО НОТ GIVE CLEAR INDICATION OF THE ORDER OF THE 


SYSTEM, INDICATING FOUR POLES WHEN THERE ARE THREE. THE 
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ZERO IN THE RIGHT HALF PLANE SHOULD BE NEGLECTED. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (с) 


SYSTEM TRANSFER FUNCTION 


POLES 1 КЕАІ, IMAGINARY 
1 -1 . 1300231 0.0 J 
2 -10.158884 0.0 J 
3 -1181775 - 536 285 012 J 
4 - 11088 7 25 +53. 205012 J 
ZEROES | REAL IMAGINARY 


GAIN CONSTANT = 252014.00 


SYSTEM STATE VARIABLE (PHASE FORM) 


1 ° 


A VECTOR 
1 34070.104 
2 33763. 954 
3 3235. 2996 
4 33.964457 
B VECTOR 
4 252014.00 
C VECTOR 
1 1.0000000 
REMARKS: 


ASK FOR FOUR POLES NO ZEROES. 


RESULTS DO NOT GIVE CLEAR INDICATION OF ORDER OF THE 


SYSTEM INDICATING FOUR POLES WHEN THERE ARE THREE. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (d) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 А > 0.0 
2 =10.190%50 0:0 
3 0.16021282 29. 91990 
4 0. 165021252 59.917890 
ZEROES REAL IMAGINARY 
1 j; 95959 73 — 02.102372 
2 3225557973 SOTO 


GAIN CONSTANT = 332.37549 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 40965.288 
2 40602.813 
3 3597.9654 
4 10.990004 
B VECTOR 
4 332. 37549 
C VECTOR 
1 ЭЛОС БЛӘ 
2 -7.1115946 
3 1.00000000 


Se 
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REMARKS: 
ASK FOR FOUR POLES TWO ZEROES 


RESULTS INDICATE TWO POLES WHEN THERE ARE THREE. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (е) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0459484 0.0 J 
2 = F0 257872 0.0 J 
3 2.2244020 0.0 J 
ZEROES REAL IMAGINARY 
1 2242035545 0.0 J 


GAIN CONSTANT = 83.810272 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 23.852104 
2 SS 8555922 
3 13.522222 
B VECTOR 
3 807S 10272 
C VECTOR 
1 24,208 3995 
2 1.0000000 
REMARKS: 


ASK FOR THREE POLES ONE ZERO. 
RESULTS INDICATE TWO POLES WHEN THERE ARE THREE. 


POLE ATE-2.2 IS CANCELLED BY ZERO AT -2.22. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (Е) 


ж” 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 =1. 0934113 0.0 Ј 
2 =10. 207312 0.0 J 


ZEROES REAL IMAGINARY 
GAIN CONSTANT = 90.3204350 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

1 11.232745 

2 11.300723 
B VECTOR 

2 90.320435 
C VECTOR 

1 1.0000000 

REMARKS: 


ASK FOR TWO POLES NO ZEROES. 
RESULTS INDICATE TWO POLES NO ZEROES WHEN THERE ARE 
ВЕ. METHOD NEGLECTED POLE AT -100.0, OTHERWISE IDENTI- 


FOTCATION VERY GOOD. 
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EXAMPLE 3 


SYSTEM TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 ле = J 
2 MS 51.5 y 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 2950.0 


TESTS FOR EXAMPLE 3 


RUN NO. OF POLES NO. OF RESULTS RESULTS 


a -2 0 Good 


D.C. POSITIONING SERVO FREQUENCY RESPONSE IN FIG. 2 
INPUT FUNCTION STEP-WISE INCREASING SIMILAR TO THE ONE 
SHOWN IN FIG. 1, BUT WITH STEP AMPLITUDE ABOUT 0.5 VOLTS. 


INITIAL CONDITIONS EQUAL ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 3 саба аа 
2 -5.3585933 5062771299 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 3332,2095 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

I 3251.6948 

2 60»; 727187 
В МЕСТОВ 

2 333222095 
C VECTOR 

1 1.0000000 

REMARKS: 


ASK FOR TWO POLES NO ZEROES. 


GOOD RESULTS. 
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EXAMPLE 4 


SYSTEM ТО BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -4.9 0.0 Ј 
2 c 0.0 J 


ZEROES REAL IMAGINARY 


GAIN CONSTANT 15090 


Ц 


TESTS FOR EXAMPLE 4 


RUN NOS OF POLES NO. OF ZEROES RESULTS 


a 2 0 Good 


DC. “ESSITIONED SERVO. 
FREQUENCY RESPONSE IN FIG. 3. 
INPUT FUNCTION THE SAME AS IN EXAMPLE 2. 


INITIAL CONDITIONS EQUAL ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -5.2129680 0.0 
2 -36.525846 0.0 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 178.37120 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

1 190.40807 

2 41.738815 
B VECTOR 

2 176.37120 
СУМЕСТОВ 

1 1.0000000 

REMARKS: 


ASK FOR TWO POLES NO ZEROES. 


GOOD RESULTS. 


60 








EXAMPLE 5 


SYSTEM TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0 0.0 Ј 
2 -10.0 0.0 J 
3 < OREO 0.0 J 
ZEROES REAL IMAGINARY 
1 ^ 5.0 0.0 Ј 


GAIN CONSTANT = 500.0 


TESTS FOR EXAMPLE 5 


RUN NO. OF POLES NO. OF ZEROES RESULTS 
a 3 1 Good 
b 2 0 See Table of Data 


INPUT FUNCTION STEP-WISE STRICTLY INCREASING SIMILAR TO 
THE ONE SHOWN IN FIG. 1. 


INITIAL CONDITIONS EQUAL ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTBM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
l -0.89534766 0.0 
2 -8.3914174 0:0 
3 -88.600646 0:0 
ZEROES | REAL IMAGINARY 
1 -3.8670622 0.0 


GAIN CONSTANT = 706.88770 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 “6652677256: 
2 830.32662 
3 97.887411 
B VECTOR 
3 706.88770 
C VECTOR 
1 3.8670622 
2 1.0000000 
REMARKS: 


ASK FOR THREE POLES ONE ZERO. 
RESULTS SHOW POOR IDENTIFICATION OF THE HIGHER 


BEBOUENCY POEE,. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -3.2441577 0.0 J 
2 - 1088917738 0.0 J 


ZEROES REAL IMAGINARY 
GAIN CONSTANT = 763.43213 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

1 3512729309 

2 111.66154 
B VECTOR 

2 763243243 
С МЕСТОВ 

1 1.0000000 

REMARKS: 


ASK FOR TWO POLES NO ZEROES. 
RESULTS CAN BE CONSIDERED AS A LOWER ORDER MATHEMATICAL 


MODEL. 
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ЕХАМРЉЕ 6 


SYSTEM TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 = 0.0 J 
2 гід Zoo 7 
3 о 0.0 5 
4 -100.0 0.0 
ZEROES REAL IMAGINARY 
GAIN CONSTANT = 660000 
TESTS FOR EXAMPLE 6 
RUN NO. OF POLES NO. OF ZEROES RESULTS 
Bad 
Good 


See Table of Data 
See Table of Data 
See Table of Data 
See Table of Data 
See Table of Data 
See Table of Data 


уат о оо o mp 
Q (Q » U о ы D 
O P P N N N O O 


INPUT FUNCTION FOR RUN (a) STEP-WISE STRICTLY INCREASING 
SHOWN IN FIG. 4, 


INPUT FUNCTION FOR RUNS (b), (c), (d), (e), (f), (g), amd 


64 








(h) STEP-WISE INCREASING-DECREASING SHOWN IN FIG. 5. 


INITIAL CONCLUSION EQUAL ZERO. 
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INPUT 
VOLTS 


80 


O TIME - MSEC 400” 


FIRST INPUT FUNCTION USED IN EXAMPLE NO. 6 


OUTPUT 
VOLTS 


60 


2 TIME - MSEC 400 
СЫ CO To Ian 220% Ti PUT: 


FIGURE 4 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.5012900 -4.3496191 Ј 
2 =1. SOLS 900 4.3496191 Ј 
3 “ =571821332 -26.089080 J 
4 -59 492137337 26.089080 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT - 62546.430 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
l ` 15047.75 
2 2366. 2077 
3 764.79050 
4 13.967046 
B VECTOR 
4 62546.430 
C VECTOR 
1 1.0000000 
REMARKS: 


INPUT FUNCTION - STEP-WISE INCREASING FUNCTION SHOWS IN 
FIG. 4. 


ASK FOR FOUR POLES NO ZEROES. 
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RESULTS NOT GOOD. 


Nene ІНЕ POSTTIVE SECOND DERIVATIVE OE. THE OUTPUT AT 


THE END OF THE SAMPLING INTERVAL. 
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р TIME - MSEC 400 


SECOHD INPUT FUNCTION USED LI BAAKPLS NO. 6 


OUTPUT 
VOLTS 
80 
P TIME - MSEC 400 


SLOTPUT COROT TNT TO Tas ABOV 11 PUT 


FIGURE 5 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 =2 0192066 OO J 
2 -10.764637 0.0 J 
3 =3 3596973 0.0 J 
4 =105. 65671 0.0 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 612357.94 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 77080.941' 
2 48396.117 
3 55512 57829 
4 152.0620 5 
B VECTOR 
4 6125357094 
C VECTOR 
1 1.0000000 
REMARKS: 


INPUT FUNCTION STEP-WISE INCREASING-DECREASING SHOWN 
ШІ EIG, 5. 


ASK FOR FOUR POLES NO ZEROES. 
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RESULTS VERY GOOD. 
NOTE THE NEGATIVE SECOND DERIVATIVE DURING SOME TIME 


WITHIN THE SAMPLING INTERVAL. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 


SYSTEM TRANSFER FUNCTION 


POLES 


1 


ZEROES 
1 


2 


GAIN CONSTANT 


REAL 
-2.0375228 
-10.649015 
SN. 2 58558 
=101 98991 

6.7412854 


6.7412854 


REAL 
6770950645 


6.2635645 


= 603327.00 


IMAGINARY 


=96% 72165 


9207372163 


IMAGINARY 


2909 9] NI» 


-96.069197 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 


1 


B VECTOR 


6 


0.7075898D 
0.44024818D 
48381252.0 
1366621.4 
12569.020 


135.45544 


603327.00 


72 


09 


09 


q 


с 


ч q GQ с 





СЕС ТОК 


1 9 2108 моею 

2 O. S2129 

3 1.0000000 
REMARKS: 


ASK SIX POLES TWO ZEROES. 
RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 


STEM, SINCE COMPLEX POLES ARE CANCELLED BY COMPLEX ERRORS. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (а) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 

1 -2.0305074 0.0 

2 -10.677930 0.0 

3 234.303456 0.0 

4 | мас 5 0.0 
ZEROES REAL IMAGINARY 

1 216.60587 0.0 

2 751088224 0.0 
GAIN CONSTANT = -1.4570827 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 71464.878 
2 44715.215 
3 4974.7267 
4 143.08956 
B VECTOR 
4 -1.4570827 
C VECTOR 
1 -388400.95 
2 =2 39322636 
3 1.0000000 
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REMARKS : 

ASK FOR FOUR POLES TWO ZEROES. 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
SYSTEM, SINCE THE TWO EXCESS ZEROES APPEAR ONE IN THE 


P. AND THE OTHER IN THE L.H.P. BUT CAN BE NEGLECTED. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (e) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -2.7344050 0.0 
2 -8.8026617 0.0 J 
3 -43.843968 0.0 E 
4 -68.587703 | 0.0 3 
5 8.0082667 0.0 E 

ZEROES REAL IMAGINARY 
Й 7.9267041 0.0 Ј 
2 -622.31546 0.0 J 


GAIN CONSTANT = 777.31934 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 -579657%55 
2 -227126.83 
3 2957. 3550 
4 3335.5841 
5 115.96047 
B VECTOR 
5 E33 4 
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REMARKS: 

ASK FOR FIVE POLES TWO ZEROES 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
Eu BM, SINCE THE EXCESS POLE AT 8.0 IS CANCELLED BY THE 
memo AT +7.9 AND THE OTHER EXCESS ZERO AT - 622.3 CAN BE 


NEGLECTED. 
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IDENTIFICATION OF UNKNOWN SYSTEM ~ Run (f) 


SYSTEM TRANSFER FUNCTIONS 


POLES REAL IMAGINARY 
1 -2.0267057 0.0 J 
2 -10.704298 0.0 J 
3 A 0.0 J 
4 -101.07630 0.0 J 
ZEROES REAL IMAGINARY 
1 -4045.2488 0.0 n 


GAIN CONSTANT = 146.26987 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
l 74636.682 
2 46730.407 
3 510291 7 
4 147.84453 
B VECTOR 
4 146.26987 
C VECTOR 
1 4045.2488 
2 1.0000000 
REMARKS: 


ASK FOR FOUR POLES ONE ZERO 


RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
78 





SYSTEM, SINCE THE EXCESS ZERO ZT - 4045.2 САМ ВЕ NEGLECTED. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (а) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
ЦІ -2.1810404 0.0 
2 -10.145443 0 0 
3 =35+910038 0.0 
3 . 96.915087 0.0 
5 16.571798 0.0 
ZEROES REAL IMAGINARY 
1 16.789359 0.0 


GAIN CONSTANT = 580130.69 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR , 
1 =127061709.1 
2 -682609.54 
3 =39334.691 
4 2752-1957 
5 128.57981 
B VECTOR 
5 580130769 
C VECTOR 
1 = 1677899359 
2 1.0000000 
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REMARKS: 
ASK FOR FIVE POLES ONE ZERO. 


RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
SYSTEM, SINCE EXCESS POLE AT 16.5 IS CANCELLED BY EXCESS 


ZERO AT 16.7. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (h) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0271020 0.0 J 
2 =-15.051165 -9.0280041 Ј 
3 —-15: 05165 9.0280041 Ј 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 4051.1528 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 316.39101 
2 |». 338.96060 
3 31.129433 
B VECTOR 
3 4051.1528 
C VECTOR 
T 1.0000000 
REMARKS: 


ASK FOR THREE PORES NO” ZEROES. 


RESULTS CAN BE CONSIDERED AS A LOWER ORDER MATHEMATICAL 


MODEL. 
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EXAMPLE 7 
ОРО TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 


1 -25.0 0.0 J 

2 za 0.0 E 

3 ESO SENS J 

4 21507 185.3 7 
2ЕКОЕ5 REAL IMAGINARY 

1 =] 0.0 E 

2 00 0.0 B 


TESTS FOR EXAMPLE 7 


RUN NO. OF POLES NO. OF ZEROES RESULTS 
a 6 3 See Table of Data 
b 4 2 Good 
C 4 3 See Table of Data 


MEUT FUNCTION SIMILAR TO THE ONE SHOWN IN FIG. 5. 


INTTIAL CONDITIONS EQUAL TO ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTEM TRANSFER FUNCTION 


ча ааа q 


POLES REAL IMAGINARY 
1 -26.773050 0.0 
2 -72.749634 0.0 
3 -101.14084 173.81324 
4 -101.14084 -173.81324 
5 960.71376 0.0 
6 52.080470 0.0 
ZEROES REAL IMAGINARY 
1 50.693422 0.0 
2 -3.0909162 0.0 
3 -9,5071748 0.0 
GAIN CONSTANT = -0.33196749E 09 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 0.39410705D 13 
2 0.14131408D 12 
3 -0.12683604D 10 
4 -43800405. 
5 -193111.44 
6 =/10 93986 

B VECTOR 
6 =0: 33196749E 09 
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C VECTOR 


1 ЕЮ. 6708 

2 -609.25446 

E 2392095331 

4 1.0000000 
REMARKS: 


ASK FOR SIX POLES THRE ZEROES. 
RESULTS INDICATE FOUR POLES TWO ZEROES, SINCE POLE 
AT 52.0 IS CANCELLED BY ZERO AT 50.7 AND POLE AT 960.7 


CAN BE NEGLECTED. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 


SYSTEM TRANSFER FUNCTION 


POLES 


1 


ZEROES 
1 
2 


GAIN CONSTANT 


REAL 
-26.992965 
АЛ Ов 
-126.,52085 


-T26.52905 


REAL 
-3.1577742 


-9.376548. 


= 444891.63 


IMAGINARY 
0.0 
0.0 
савана вын 


1522989579 


IMAGINARY 


SYSTEM STATE VARIABLES (PHASE FORM) 


А МЕСТОК 
1 
2 
3 
4 

B VECTOR 
4 


C VECTOR 


0.10367784D 
5706927.6 
78054.592 


357.64263 


444891.63 


291700902 2 


12.334322 


71.0000000 
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REMARKS: 


ASK БОК FOUR POLES TWO ZEROES. 


GOOD RESULTS. 
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IDENTIFICATION OF UNENOWN SYSTEM - Run (с) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 27 3858007 0:40 
2 -75.010363 02% 
3 -119, 39880 =181 29378 
4 = 119. 37550 181.29378 
ZEROES REAL IMAGINARY 
J =3 71678969 0.0 
2 -9,4693811 0.0 
3 -5036.4256 0.0 


GAIN CONSTANT = 82.045258 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 96790162.0 
2 531 800.5 
3 73620.633 
4 341.15318 
B VECTOR 
4 82.045258 
C VECTOR 
1 14959912 
2 63519.731 
3 ` 5049.0317 
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3 1.0000000 


REMARKS: 
ASK FOR FOUR POLES THREE ZEROES. 
RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 


SYSTEM, SINCE ZERO AT -5036.4 CAN BE NEGLECTED. 
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EXAMPLE 8 
Seel TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0 0.0 J 
2 -10.0 0.0 J 
3 25.0 0.0 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 2500 


TESTS FOR EXAMPLE 8 


RUN NOW Or POLES NO. OF ZEROES RESULTS 
a 3 0 Good 
b 3 0 “Good 
С 3 0 Good 


INPUT FUNCTION SIMILAR TO THE ONE SHOWN IN FIG. 5. 
INITIAL CONDITIONS DIFFERENT FROM ZERO FOR THE THREE 


DIFFERENT RUNS. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
iL -0.97837856 0.0 
2 S10 MOOT 0.0 
3 -24.694433 0.0 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 2378.9880 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 2655209877 
2 305% 93838 
3 36.649709 
В МЕСТОК 
3 2378 298180 
C VECTOR 
1 1.0000000 
REMARKS : 


ASK FOR THREE POLES NO ZEROES. 


GOOD RESULTS WIPRA” I.C. DIPPERENT FROM" ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
j ] -0.96861216 0.0 
2 -10.956655 0.0 
3 224.730535 0.0 

ZEROES REAL IMAGINARY 


GAIN CONSTANT = 2379.6726 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 263.01082 
2 306. 15109 
3 36.707801 
B VECTOR 
3 2379.6726 
C VECTOR 
1 1.0000000 
REMARKS: 


AORN FOR THREE PORES NO ZEROES. 


COSPSREUSUETS WITH Т.С, DIFFERENT FROM ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 a24. 839515 0.0 J 
2 - 4950069923 0.0 J 
3 - 10297 3508 0.0 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 2386.2974 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 261.64609 
2 s 306295352 
3 Ji 772923 
B VECTOR 
3 2386.2974 
C VECTOR 
1 1.0000000 
REMARKS: 


ASE BOR THREE POLES NOI ZEROES. 


GOOD RESUBTS WITH I.C. DIFFERENT FROM ZERO. 
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EXAMPLE 9 


DeNOWEM"TO BE IDENTIFIED 


SAME SYSTEM EXAMPLE 8. 


TESTS FOR EXAMPLE 9 


RUN NO. OF POTES NO. OF ZEROES RESULTS 

a 3 0 Bad 

b 3 0 Almost Good 
e 3 0 Good 

а 3 0 Good 


INPUT P'UNCTION STEPSWISE INCREASING-DECREASING SIMILAR 
TO THE ONE SHOWN IN FIG. 5. 

INITIAL CONDITIONS EQUAL ZERO. 

LOW FREQUENCY NOISE - FROM A FEW CYCLES TO ABOUT 300. 
CARES FER SECOND = WAS ADDED AT THE INPUT AND OUTPUT. 


БЕТТІ R.M.S. VALUES OF NOISE WERE CONSIDERED., 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -0.257856 14D-01 0.0 J 
2 -16.598045 -7.2418627 J 
3 160598045 7. 2А 086207 Ј 


ZEROES REAL IMAGINARY 


GAIN CONSTANT 2294.0054 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 8.4561459 
2 3218. 2856 / 
3 33 , 2287 
B VECTOR 
8 2294.0054 
С VEGEOR 
1 1.0000000 
REMARKS: 


ASK FOR THREE POLES NO ZEROES. 
Ban 5. 2NVZEUEFSORZTEEZNOTSE,1.0.VOLTS. 


RESULTS NOT GOOD. 
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IDENSIEICATION OF UNKNOWN SYSTEM - RUN (b) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 —2 900/0157 0.0 
2 = ЩО КЮ З 0.0 
5 -0.72464570 0.0 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 2364.4692 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 206.93154 
2  311.29595 
3 36.236600 
B VECTOR 
3 2364.4692 
C VECTOR 
i 1.0000000 
REMARKS: 


Row FOR THREE РОШЕ NO ZEROBS. 
HOMES. VALUE OPUPHETNIOTSE . 7/35 "VOLTS. 


RESULTS ALMOST GOOD. 
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TDENBIEITCATION OP UNKNOWN SYSTEM - Run (c) 


SYSTEM TRANSEER FUNCTION 


POLES REAL IMAGINARY 
1 -1; 1475703 0.0 
2 -10.808421 0.0 
3 —29 4175558 ОШО 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 2289.9822 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 290. 17727 
2 293.07699 ` 
3 35.431549 
B VECTOR 
3 2289.9822 
C VECTOR 
1 1.0000000 
REMARKS: 


ASK FOR THREE POLES NO ZEROES. 
RS. VALUE ОР TAE NOISE 0.5 VOLTS. 


RESULTS GOOD. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (d) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0495050 0.0 J 
2 -10.680594 0.0 J 
3 -25702 3091 00 Ј 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 2381.8215 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 280.49226 
2 ` 304.73267 
3 36. 7999/0 
B VECTOR 
3 23915 215 
CEV ECTOR 
18 1.0000000 
REMARKS: 


ASK FOR THREE PULES NO ZEROES. 
ENES: VALUE OF THE NOISE .25 VOLTS. 


RESULTS VERY GOOD. 
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EXAMPLE 10 


SYSTEM TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0 0.0 7 
2 -10.0 0.0 7 
3 -20.0 0.0 Ј 
4 -33.0 0.0 Ј 
5 -100.0 0.0 Ј 
ZEROES REAL IMAGINARY 
1 -5.0 0.0 


GAIN CONSTANT = 2640000.0 


TESTS FOR EXAMPLE 10 


RUN NO. OF “ROLES NO. OF ZEROES RESULT 

a 5 L - See Table of Data 
b 3 0 See Table of Data 
С 2 0 See Table of Data 


INPUT FUNCTION SIMILAR TO THE ONE SHOWN IN FIG. 5. 
INITIAL CONDITIONS EQUAL ZERO, 
SEE FIG. 6 COMPARING THE RESPONSES OF THEORETICAL 


AND IDENTIFIED SYSTENS TO A UNIT STEP INPUT. 








IDENTIFICATION OF UNKNOWN SYSTEM - Run (а) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -0.94505946 0.0 
2 -10.807637 0.0 
3 -32.034223 0.0 
4 =62 3699979 0.0 
5 СОКО 0.0 
ZEROES REAL IMAGINARY 
1 -4.9340575 0.0 


GAIN CONSTANT = 2741280.0 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 248.538.7 
2 2992 97 L5 
3 403446.21 
4 : 16013.272 
D 227.45465 
B VECTOR 
5 2741280.0 
C VECTOR 
1 4.9340575 
2 1.0000000 
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REMARKS : 
ASK FOR FIVE POLES ONE ZERO. 
RESULTS NOT VERY GOOD FOR THE HIGHER FREQUENCY POLES. 
SEE FIG. 6 COMPARING THE RESPONSES OF THEORETICAL AND 


BENT IF MED SYSTEMS ТО А UNIT STEP INPUT. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -2.6408186 0.0 J 
2 -24.059549 0.0 Ј 
3 -24.059549 0.0 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT - 10410.082 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 3331.3020 
2 1388.5395 ` 
3 50.759916 
B VECTOR 
3 10410.082 
C VECTOR 
1 1.0000000 
REMARKS: 


ASK FOR THREE POLES NO ZEROES. 

RESULTS CAN BE CONSIDERED A LOWER ORDER MATHEMATICAL 
MODEL. 

SEE FIG. 6 COMPARING THE RESPONSES OF THEORETICAL AND 


DIT TED SYSTEMS TO A UNIT STEP INPUT. 
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TBEITTIEICADTION OF UNKNOWN SYSTEM - Run (c) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 =576608 772 0.0 J 
2 -1 3. 680660 0.0 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 182.89790 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

1 77.444535 

2 1959990537 
B VECTOR 

2 182.89790 
C VECTOR 

1 1. 0000000 

REMARKS: 


ASKIFOR TWO POLES NO. ZEROES. 

RESULTS CAN BE CONSIDERED A LOWER ORDER MATHEMATICAL 
MODEL. 

SEE FIG. 6 COMPARING RESPONSES OF THEORETICAL AND 


ТОРЛЕ О SYSTEMS TO A UNIT STEP INPUT. 
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UNTITS=OUTPEUT 





TIME-SECONDS 


UNIT STEP RESPONSE OF: 


(ал меге to be identified 

(б) Fifth order identified model 
(c) Third order identified model 
(d) Second order identified model 


FIGURE 6 
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EXAMPLE 11 


SYSTEM TO BE IDENTIFIED 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0 0.0 J 
2 -10.0 0.0 J 
3 -33.0 0.0 J 
ZEROES REAL IMAGINARY 
1 -5.0 0.0 J 
GAIN CONSTANT = 52500.0 
TESTS FOR EXAMPLE 11 
RUN NO. OF POLES NO. OF ZEROES RESULTS 
a 4 1 Good 
b 3 0 See Table of Data 
C 2 0 See Table of Data 


INPUT EUNCTION IS SIMTLCAR TO THE ONE SHOWN IN FIG. 5. 
INITIAL CONDITIONS EQUAL ZERO. 


SEE FIG. 7 & 8 COMPARING THE RESPONSES OF THEORETICAL 


ZIDERDEIIDTEYINGISYSTEMS TO ASSTEP INPUT. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 


SY STEM LRANSEER FUNCTION 


POLES REAL IMAGINARY 
1 -1.0182594 0.0 Ј 
2 -1 0.830145 0.0 J 
3 -35.841486 0.0 J 
4 =1 08727539 0.0 J 
ZEROES REAL IMAGINARY 


1 =5 70491588 0.0 J 
GAIN CONSTANT = 64711.781 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 428677999 
2 47640.814 
3 5600.2860 
= 159.9 3228 
B VECTOR 
4 64711.781 
C VECTOR 
1 5.0491588 
2 1.0000000 
REMARKS: 


ASK FOR FOUR POLES ONE ZERO. 


RESULTS GOOD. 
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SEE FIG. 7 & 8 COMPARING THE RESPONSES OF THEORETICAL 


AND IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 


* 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY | 
1 -31.063053 -36 5868 Ј 
2 -31.063053 36.791868 J 
3 -2.9265019 0.0 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 28337.477 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 
1 678 5.2550 
2 · 2500.3670 © 
3 65.052607 
B VECTOR 
3 28337.477 
C VECTOR 
1 1.0000000 
REMARKS: 


ASK FOR THREE POLES NO ZEROES. 

RESULTS CAN BE CONSIDERED A LOWER ORDER MATHEMATICAL 
MODEL. 

SEE FIG. 7 & 8 COMPARING THE RESPONSES OF THEORETICAL 


ND IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 


SYSTEM TRANSFER FUNCTION 


POLES REAL IMAGINARY 
1 =353. 438 / 3 0.0 J 
2 -3607703934 0.0 J 
ZEROES REAL IMAGINARY 


GAIN CONSTANT = 537.30396 


SYSTEM STATE VARIABLES (PHASE FORM) 


A VECTOR 

il 125.94514 

2 40.134951 
B VECTOR 

2 537. 50396 
C VECTOR 

1 1.0000000 

REMARKS: 


ASK FOR TWO POLES NO ZEROES. 
RESULTS CAN BE CONSIDERED A LOWER ORDER MATHEMATICAL 


MODEL. 
SEE IGS. 7 &6 8 COMPARING THE “RESPONSES OF THEORETICAL 


AND IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 
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IV. CONCLUSIONS 


The method of multiple integrations is a practical method 
for identifying lumped parameter, linear, time invariant 
systems. It works well in the absence of noise up to a fourth 
order system if the input function is carefully selected. 

For higher order systems the identification is in general not 
good, but an approximation of the system can be obtained. 

If desired, the method can also be used to obtain a reduced 
order mathematical model of the system, after the 

raentificat ten. 

The accuracy of an identification depends very strongly 
ра Епе nose present. Small R.M.S. values of noise are 
enough to make the identification bad. Noise with a R.M.S. 
value of about 0.5 volts and below do not practically affect 
the identification, when inputs like the ones ued in this 
research are present, independently of where the noise is 
present - input, output or both -. Accuracy also depends 
very much on the type of input function used to drive the 
system. 

The computational requirements of the method are not 
excessive, however, an accurate A/D converter system is very 
desirable and the better this operation is done, the more 
accurate will the identification be, since the accuracy of 
an identification is usually comparable to the number of 


exact significant digits in the input-output data. 
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The limitation of the method to low order systems is 
primarily due to the algorithm used to solve the over-deter- 
mined set of linear equations. As better algorithms become 
available, it will be possible to identify higher "и 
Systems. 

The results obtained when reduction of the system order 
was desired were not very good at first sight, but when the 
steady state gain of the lower order models was adjusted, 
making it equal to the one of the non-reduced identified 
model, the approximation can be considered good. However, 
when a pair of dominant complex roots is present in the 


system a better reduced order model should be expected. 
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APPENDIX I 


EN OCTET GO THE INPUT = OUTPUT RECORD 

The identification technique described in Chapter II is 
based on knowledge of the input and output record and since 
the IBM 360 Digital Computer was used to perform the inte- 
grations of equation (5), the record should be available in 
a digital form. Thus, and A/D (analog to digital) operation 
should be performed. This A/D operation was accomplished 
using the Hybrid Computer System (XDS 9300/Ci 5000) located 
in the Electrical Engineering Computer Laboratory. The 
output of this computer is a seven track magnetic tape in 
which the input and output functions are recorded in digital 
form. Since the XDS 9300/Ci 5000 system m used on a self- 
service basis, a brief discussion of what was done to obtain 
the record, is given here, so that the same procedure could 
be used in the future whenever an A/D conversion is needed - 
it should be pointed out here than the computer center staff 
was always predisposed to assist in the operation of the 
system. 

The analog part of the system, the Ci 5000, has two patch- 
boards which are the heart of the system. The logic patch- 
board was utilized to set the sampling rate, at which the 
analog signal was sampled - 10KHZ in this particular case. 

This sampling rate is the maximum that can be achieved 
with this computer, and can be changed (decreased) by 


altering a counter-like number display located in the 
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frontal panel of the computer - FP00 and FP04. Also in this 
patchboard the base frequency - 1KHZ - used to determine the 
moments at which the input function is changed step-wise, is 
set. The same counter-like number display can be used to 
change this frequency - FF10 and 14. In the digital patch- 
board, the wiring that controls the digital part of the A/D 
operation is also set. 

The analog patchboard was used to introduce the proper 
isolation between the different components of the different 
systems tested and to provide the desired amplification 
factor. This patchboard was also used in the generation of 
the input Signal. This signal was obtained in the following 
way. An oscillator generates a Sine wave in which amplitude 
and frequency can be adjusted by two potentiometers - P043 
and P417 -; the sine-wave is amplified and passed through a 
zero order hold circuit which output is a step-wise wave form, 
that was used as the input function. The sine wave was also 
used to control the beginning of the sampling. The wiring of 
both patchboards was set by the computer center staff and 
those patchboards - numbered 21 and 21 - are reserved for 
the A/D operation. 

The input and output Signals are fed from the analog 
part of the computer system to the digital part. The 
Digital part interfaced with the analog is a Xerox Data 
System - XDS9300. 

Two tapes drive units are available as input-output 


devices. One printer is used for program listings and 
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21. With the sense switches at the XDS9300 console off, 
depress CLEAR and CLEAR FLAGS buttons simultaneously. 

22. Mount the magnetic tape following the instructions 
available on the inside of the protective doors. 

23. Select density 556 bits per inch with the correspond- 
ing knob on the face of the tape driver. 

24. Read in the A/D program, described in Chapter III 
and listed in Chapter IV, following the instructions available 
on the console of the XDS9300. 

25. Enter, through the teletype, the parameters to be 
used for the A/D program. This program has seven options or 
instructions ака entered also through the teletype by 
a single digit number and depressing the carrlage return 
button. If a further response is required the computer would 
indicate the format. The options are: 

1 = Desire to enter new parameters. 
2 = Start digitizing the analog input signals 
(normally actuated by manual switch operation on Ci 5000 


control panel; (р51 up and key COMPUTE depressed). 


3 = Write an END OF FILE on the magnetic tape. 

4 = Rewind the tape to the load point (beginning). 
5 = Эр the following number ofsMND OFWBELE marks. 
6 = Print the following amount of digitized data. 

7 = Actuate the A/D program for the next Single 


block of data encountered and feed the signal to the strip 
chame reeord@er (Witich mW@st “se running prior te this option 


selection connected to the analog patchboard on the Ci 5000. 
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merest of the Counters should be set to Zero. 

7. Depress the master clock control, RUN. 

Ө. SET 00 (Delay flip-flop) to .l ms. and adjust it, 
until its corresponding light on the front panel of the 
Ci 5000 comes on. 

9. Set DFO6 in .1ms., DFO7 in 1 sec. and DF10 in 1 ms. 

10. Set all the DS (digital switches) to the center 
position. 

I1. Torn on the TRC. 

12. Set switches Xl to the left and X2 to the right. 

13. Select amplifier 042 in the Ci 5000 console, adjust- 
ing P417, DF06 and F10, F14 to obtain the desired input 
function which сап be observed in the TRC when the DS1 
switch is set in the up position and the button COMPUTE 
is depressed in the Ci 5000 console. 

AS SCA Ree. СТЕК switch on XDS9300 to EXIT. 

15. If the XDS9300 is not energized, depress ‘simultaneously 
RESET and POWER buttons. 

16. If the teletype is not ON, energize it by turning the 
knurled knob Со ON. 

17. If hte Line Printer is not ON, depress the POWER 
switch - indicator on the face of the printer. About one 
minute later the lower half of the Power button/lights up 
indicating the printer is energized. 

18. Depress READY SWITCH INDICATOR. 

19. Depress the POWER ON button on the card reader. 


20. Depress IDLE button on the console of the XDS9300. 
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outputs. A card reader is associated with the digital com- 
puter as the avenue of input for compiling programs. A 
teletype unit is used for parameter.inputs and instructions 
to the digital computer program. A control console displays 
information of the computer and has the corresponding control 
switches. A recorder is also available with which the 
records of the input and output, after the A/D operation is 
finished, can be recorded in a strip chart. 

After the preceding general description of the system, 
the following procedure can be followed to take a record. 

1. Set the analog and digital patchboards in the 
Ci 5000. 

2. Connect the output of the amplifier 042 to the 
input of the system. 

З. “CGemmect the output of the system to the input of 
the amplifier 036. 

4. If the Ci 5000 is not energized, depress on switch 
indicator at lower left corner of display panel of Ci 5000. 
Depress KEY BOARD and POTSET buttons on the display panel. 

5. Enter the maximum value of the sine wave through 
the POTSET, POT and SERVO keys. 

6. Set the counters FF's in the following positions: 

ЕЕ00 іп 0, ЕЕ04 in 9, for 10KHZ sampling rate. 

FF10 in 2, F14 in 4, for 40HZ rate of changing, 
step-wise, the input function. 

FF30 in O, if only one record per run is desired; 


increasing FF30 the number of records/run can be adjusted. 
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The parameters to be entered in point 25 are: 


NSAMP = Total number of samples to be taken from 
each channel - 4002 in this case. 
NCHAN = Number of channels to be digitized - 2 in 


this case. 

| NREC = Number of records to be digitized per run; 
this number should be equal to the one set in the counter FF30 
(Ci 5000) plus 1; in this particular problem NREC = 1. 

ITAPE = The identifying number of the tape drive in 
use. This number was set when the tape was mounted. 

NDEL = Controls the delay between the digital and 
analog signal in option 7. 

After entering the parameters the * character should be 
typed on the teletype followed by RETURN. At this point 
everything is ready to take a record; by typing, 2 - RETURN, 
on the teletype, and setting switch DS1 - UP and pressing 
COMPUTE on the Ci 5000. After the record was taken, write 
- 3 - (END OF FILE). If it is desired, the record can be 


recorded in a strip chart recorder. 


В. SEVEN TO NINE TRACK TAPE CONVERSION 

In the magnetic tape generated in point A, the input- 
output record is written in a seven track tape, but since 
this tape is going to be used with the IBM 360 Digital 
Computer and this computer uses NINE TRACK TAPES, a 
conversion of the seven track tape into a nine track tape 


is required. 
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This operation is done by the IBM 360, using the program 
CONVERT, which is described in the Technical Note No. 0211-08 
second edition - PROCEDURES FOR CONVERTING 7-TRACK MAGNETIC 
TAPES TO 9-TRACK MAGNETIC TAPES by Sharon D. Raney [21]. 

This description is available at the Computer Center in 
Ingersol Hall. A listing of the CONVERT program is given 


in Chapter IV. 
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CONVERT 
THIS PROGRAM CONVERTS A SEVEN TRACK TAPE TO A NINE TRACK TAPE 
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